Effect of a neutron-star crust on the r-mode instability 
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The presence of a viscous boundary layer under the solid crust of a neutron star dramatically 
increases the viscous damping rate of the fluid r-modes. We improve previous estimates of this 
damping rate by including the effect of the Coriolis force on the boundary-layer eigenfunction and 
by using more realistic neutron-star models. If the crust is assumed to be perfectly rigid, the 
gravitational radiation driven instability in the r-modes is completely suppressed in neutron stars 
colder than about 1.5 x 10 8 K. Energy generation in the boundary layer will heat the star, and will 
even melt the crust if the amplitude of the r-mode is large enough. We solve the heat equation 
explicitly (including the effects of thermal conduction and neutrino emission) and find that the r- 
mode amplitude needed to melt the crust is ct c ~ 5 x 10~ 3 for maximally rotating neutron stars. If 
the r-mode saturates at an amplitude larger than a c , the heat generated is sufficient to maintain 
the outer layers of the star in a mixed fluid-solid state analogous to the pack ice on the fringes of 
the Arctic Ocean. We argue that in young, rapidly rotating neutron stars this effect considerably 
delays the formation of the crust. By considering the dissipation in the ice flow, we show that the 
final spin frequency of stars with r-mode amplitude of order unity is close to the value estimated 
for fluid stars without a crust. 



I. INTRODUCTION 

The r-modes (fluid oscillations governed primarily by 
the Coriolis force) have been the focus of considerable at- 
tention over the past few years (see Friedman and Lock- 
itch [1] for a review). The gravitational-radiation driven 
instability of these modes has been proposed as an ex- 
planation for the observed relatively low spin frequen- 
cies of young neutron stars, and of accreting neutron 
stars in low-mass x-ray binaries (LMXBs) as well. The 
r-mode instability may also provide a source of gravi- 
tational waves detectable by the "enhanced" LIGO and 
VIRGO interferometers, the configurations whose opera- 
tion is expected to begin in about the year 2006. In real 
neutron stars this instability can only occur when the 
gravitational-radiation driving timescale of the r-mode 
is shorter than the timescales of the various internal dis- 
sipation processes that may occur in the star. In this 
paper we re-examine and re-calculate some of the dissi- 
pation timescales associated with the crust of the neutron 
star. 

Recently Bildsten and Ushomirsky [2] made the first 
estimate of the effect of a solid crust on the r-mode in- 
stability. They found that the shear dissipation in the 
viscous boundary layer between the solid crust and the 
fluid core decreases the viscous damping timescale by 
more than a factor of 10 5 in old, accreting neutron stars 
and more than 10 7 in hot, young neutron stars. The 
viscous damping timescale in these stars is thus com- 
parable to the gravitational radiation driving timescale, 
and so Bildsten and Ushomirsky concluded that the r- 
mode instability is unlikely to play a role in old, accret- 
ing neutron stars. In hot, young neutron stars they also 



predicted that this boundary-layer damping mechanism 
severely limits the ability of the r-mode instability to re- 
duce the angular momentum of the star, and hence to 
produce detectable amounts of gravitational radiation. 

However, the debate over the relevance of the r-mode 
instability to the observed spin periods of neutron stars is 
not settled. Andersson et al. [3] corrected a minor numer- 
ical factor in the work of Bildsten and Ushomirsky [2] and 
used different neutron-star parameters to obtain a signif- 
icantly different result for the critical frequency of the 
onset of the r-mode instability. Andersson et al. [3] esti- 
mated the critical frequency (the frequency of rotation of 
the star at which the driving and damping timescales are 
equal) to be about 40% lower than the estimate of Bild- 
sten and Ushomirsky [2]. Based on this new estimate, 
and contrary to the conclusions of Bildsten and Usho- 
mirsky [2], Andersson et al. [3] inferred that the r-mode 
instability is likely to be the mechanism that limits the 
LMXB spin periods and those of other millisecond pul- 
sars as well. 

The calculations of Bildsten and Ushomirsky [2] and 
of Andersson et al. [3] depend on an extremely simple 
model of the boundary layer which neglects the Coriolis 
force, the dominant restoring force for the r-modes. Rieu- 
tord [4] , building on the work of Greenspan [5] , improved 
this model by finding the self-consistent solution to the 
linearized Navier-Stokes equations (including the Coriolis 
force terms) throughout the boundary layer. Rieutord's 
model of the boundary layer includes the correct angular 
structure and results in dissipative timescale estimates 
that are significantly shorter than those of Andersson et 
al. [3]. Coincidentally, Rieutord's estimates of the criti- 
cal spin frequencies agree rather closely with the original 
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estimates of Bildstcn and Ushomirsky [2] . These changes 
in the values of the boundary layer dissipation timescale, 
and the corresponding changes in the conclusions regard- 
ing the physical relevance of the r-mode instability for 
neutron stars with crusts, have prompted us to revisit 
this subject once again. 

Our aim is to provide a complete, careful re-derivation 
of recent results [2-4], including the effects of non- 
uniform density stellar models — an important factor ne- 
glected by Andersson et al. [3] and Rieutord [4]. We use 
a variety of stellar models to explore the sensitivity of 
the results to the poorly-known neutron-star equation of 
state. In summary (see Fig. 1), we find that the criti- 
cal frequency is even greater (by 25% to 50% depending 
on the equation of state) than that estimated by Bild- 
sten and Ushomirsky [2] and Rieutord [4] , and double to 
triple that estimated by Andersson et al. [3]. Our new 
results make it appear unlikely that the r-mode insta- 
bility is responsible for limiting the spin periods of the 
LMXBs. 

However, the interpretation of our results is somewhat 
complicated by the recent work of Levin and Ushomir- 
sky [6]. They showed that mechanical crust-core cou- 
pling can reduce the relative velocity between the crust 
and the core, thereby reducing the shear and the viscous 
dissipaton in the boundary layer. Let Av/v denote the 
difference between the velocities in the inner edge of the 
crust and outer edge of the core divided by the velocity 
of the core. In the static, rigid crust case Av/v = 1, but 
Levin and Ushomirsky [6] find that this quantity can lie 
anywhere in the range 0.05 < Av/v < 1. They find that 
the precise value of Av/v varies in a complicated manner 
with the spin frequency of the star, and is quite sensitive 
to the thickness of the crust and other model parame- 
ters. Throughout this paper we assume that the crust 
is static, i.e. that Av/v = 1, but indicate at appropriate 
points in the text how to rescale various quantities (such 
as dissipation timescales) for Av/v < 1. Our calculations 
therefore provide only an upper limit on the critical fre- 
quency for the onset of the r-mode instability in neutron 
stars with crusts. However our limits are easily adjusted 
once a more realistic value of Av/v is known. 

In young neutron stars, the r-mode instability is still 
a viable mechanism for spindown even if Av/v = 1. In 
the presence of a crust, the majority of the viscous dis- 
sipation is confined to a very thin boundary layer. If 
the r-mode amplitude is larger than some critical value, 
the resulting heating in this layer is so intense that it can 
compete with neutrino cooling and heat the crust-core in- 
terface to the melting temperature. This possibility was 
first suggested by Owen [7], who crudely estimated the 
critical amplitude to be a c ~ 10~ 3 , and the idea of crust 
re- melting was suggested again by Andersson et al. [3]. 
Here we perform a comprehensive, self-consistent analy- 
sis of this heating effect, including conductive transport 
of heat into the core and the crust, and eventual radiation 
of the excess thermal energy by neutrinos. We find that, 
for maximally rotating neutron stars, the critical dimen- 



sionless r-mode amplitude needed to heat the crust-core 
interface to the melting temperature at the equator is 
a c i=a 5 x 10~ 3 . [This value depends on the spin fre- 
quency and is somewhat larger away from the equator; 
see Eq. (42) and Fig. 4.] If the r-mode amplitude grows 
to a value exceeding a c in a hot, young neutron star, the 
crust will not form as usual and the neutron star will spin 
down to a much lower frequency than would have been 
possible were a crust present. 

What happens if the r-mode amplitude does exceed 
the critical value, a c l If the r-mode completely melts the 
crust, the boundary layer and the heating are removed 
and the outer layers of the neutron star quickly cool back 
down to melting temperature. If the crust completely 
cools and solidifies, the boundary-layer heating due to 
the r-mode quickly heats the crust back up to melting 
temperature. It is clear that, in the presence of a strong 
enough r-mode, neither a solid crust nor a pure fluid is 
possible. We imagine the situation to be similar to the 
pack ice on the Arctic Ocean. While the r-mode am- 
plitude exceeds the critical value a c in this picture, the 
outer layers will be composed of chunks of crustal "ice" 
floating in the fluid at the melting temperature of about 
10 10 K. The dissipation mechanism in this pack ice will be 
a combination of macroscopic viscosity due to collisions 
between chunks of ice and microscopic viscosity due to 
boundary layers bordering the chunks. Calculation of the 
viscosity in this situation would be very complicated but 
for the fact that the system must be maintained close 
to the melting temperature. The size of the ice chunks 
(and other variables controlling the viscosity) will adjust 
themselves so that the heat dissipated in the ice flow bal- 
ances the neutrino cooling. This allows us to estimate the 
r-mode damping timescale quantitatively in the pack ice, 
without knowing the details of this complicated process. 
We find that, for r-mode saturation amplitudes of order 
unity, the final spindown frequency is little changed from 
that of the purely fluid model of the instability considered 
over the past few years. 

The rest of this paper is organized as follows. In Sec. II 
we re-derive the velocity profile of the r-modes in the 
boundary layer, and in Sec. Ill we re-derive the energy 
dissipation in the boundary layer using techniques and 
notation that will be more familiar to the relativistic as- 
trophysics community. In Sec. IV we apply these results 
to the question of the stability of the r-modes in hot, 
young neutron stars and in older, colder, accreting ones. 
In Sec. V we derive the thermal structure of the bound- 
ary layer and find the r-mode amplitude a c necessary to 
heat the bottom of the crust to melting. In Sec. VI we 
argue that the presence of an r-mode with an amplitude 
greater than a c will in fact prevent crust formation, and 
instead lead to the pack-ice flow described above. We 
compute the effective dissipation in this flow, and con- 
sider the implications of the delayed crust formation for 
the developent of the r-mode instability. In Sec. VII we 
summarize and discuss some of the implications of our re- 
sults. In the Appendix we summarize the relevant ther- 
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modynamic properties of the neutron star matter near 
the crust-core interface. 



II. STRUCTURE OF THE BOUNDARY LAYER 

We begin by re-deriving the structure of the r-modes 
in the boundary layer near a rigid solid crust. Our anal- 
ysis improves the initial studies of Bildsten and Ushom- 
risky [2] and of Andersson et al. [3] by properly evalu- 
ating the angular structure of the boundary layer. We 
follow closely the more recent work of Rieutord [4], but 
employ a notation that is more familiar to the relativis- 
ts astrophysics community and improve his estimates by 
allowing non- uniform density stellar models. 

Let us decompose the fluid perturbation representing 
an r-mode into two parts: the cigenfunctions describing 
the mode in the zero viscosity limit, Sv a and SU = Sp/p, 
and the corrections that must be added to these when 
viscosity is present, Sv a and SU. The velocity correction 
Sv a must be chosen so that the relative velocity between 
the fluid core and the solid crust vanishes: for the case 
of a rigid crust this is equivalent to = Sv a + Sv a . Thus 
the viscous corrections to the velocity field are not small, 
at least near the crust. The correction to the hydrody- 
namic potential SU, however, will turn out to be small 
everywhere. 

The equations for the viscous corrections to the ve- 
locity field are obtained by expanding the Navier-Stokcs 
equation to first order in the amplitude of the pertur- 
bation. We assume that the equilibrium star is rigidly 
rotating with angular velocity ft, and that the tempo- 
ral and angular dependence of the mode is e l,jJt + tm v . As 
usual in boundary-layer theory, we assume that the fluid 
functions change much more rapidly in the direction per- 
pindicular to the boundary. Thus we assume that the 
angular derivatives are much smaller than radial deriva- 
tives, and neglect them. Under these assumptions, the 
equations that determine the viscous corrections to the 
fluid flow are: 

t(w + mtt)5i r - 20r sin 2 95i v = -d r 8U + ^d?Sv r , 
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where r] and p are the viscosity and density of the fluid, 
respectively. These equations assume only that the vis- 
cous corrections vary much more rapidly with r than with 



the angular coordinates. Thus we generalize the analy- 
sis of Rieutord [4]: we do not assume a priori that the 
fluid flow is incompressible, or that the equilibrium is of 
uniform density. 

The solutions to these equations depend exponentially 
on r with a characteristic length scale d: 



d z 



V 
2Qp 



(4) 



r=R c 



The radial velocity correction Sv r must vanish both at the 
boundary and deep within the fluid. Given the exponen- 
tial nature of the solutions, it follows that Sv r = every- 
where within the fluid. With this simplification, Eqs. (2)- 
(3) determine the velocity corrections 5v e and 8v v , while 
Eq. (1) determines SU in terms of Sv v . Eqs. (2) and (3) 
may be re-written in the following form: 



6»)] + ^81116*^^) = 0, (5) 



d 2 dt 



'(20 



cos I 



(8v e - i sin 9Sv v ) = 0, (6) 



where a = w + mil is the mode frequency in the rotating 
frame. 

It is straightforward now to write down the general so- 
lutions to these equations, and then to impose the bound- 
ary condition, = Sv a +Sv a , at the inner edge of the crust 
r = R c . These solutions are given by 

Sv e = -Sv e (R c )A+(r,9) - ism98v v (R c )A-(r,9), (7) 



sin 9 Si* = - sm9 8v ip (R c )A+(r, 9) + i8v e (R c )A-(r, 9), 

(8) 

where Sv a (R c ) is the standard non- viscous r-mode veloc- 
ity perturbation, evaluated at the radius of the inner edge 
of the crust r = R c . The functions A±(r, 9) are defined 

by 



A ± ( rj 0) = l e -(\/i(cos8+<r/2n) ± 1 e -C\A( cos 0-<r/2fi) ^ (g) 

where £ is the dimensionless radial parameter 

R c -r 



C = 



(10) 



and the characteristic thickness of the boundary layer, 
d, is defined in Eq. (4). For r-modes the frequency of 
the mode (as measured in the co-rotating frame of the 
fluid) that appears in the definition of A± has the value 
er/20 = l/(m + 1). For rapidly rotating neutron stars 
d w (10 8 K/T) cm, where T is the temperature in the 
boundary layer. Therefore, d/R c <C 1. 
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III. DISSIPATION IN THE BOUNDARY LAYER 

The shear of the r-mode velocity field is dominated 
by the rapid radial change in 8v a through the boundary 
layer. Thus, up to terms of order d/R c , the square of the 
shear tensor in the boundary layer is given by 

5a* ab 5a ab = \R 2 C (\d r Si e \ 2 + sin 2 | W| 2 ) . (11) 

The angular structure of Sv a in Eqs. (7)-(8) is deter- 
mined in part by the structure of the dissipation-free ve- 
locity field 6v a . For convenience we may write the non- 
dissipative r-mode velocity field as 



Sv e = -iAr" 1 " 1 sin™" 1 e *M+m V ) 



Sv v = Ar m - 1 sin m - 2 cos<9e 4M+mv) , 



(12) 



(13) 



where A is a normalization constant. It is straightforward 
then to evaluate the square of the shear tensor: 



I 412 2m 

Sa* ab 6<T ab = l -L^F(r,6), 



(14) 



where 



F(r, 0) = sin 2 " 1 - 2 [(1 - cos 0) 2 p* e -Cp + 

+ (l + cos6») 2 3 9 2 _e- Cp -], (15) 



and 



p± = \/2|cos0± l/(m+ 1)|. 



(16) 



Now integrate the energy dissipation rate due to shear 
viscosity over the fluid interior to the crust, ignoring 
terms of order d/R c : 



^ = - J 2r]5(T* ab 5a ab d 3 x = 2V2ir\A\ 2 ndR 2 c m+2 p c l„ 



(17) 

where l m is defined by 

T m = f sin 2 ™" 1 0(1 + cos0) 2 y/\cos0-l/{m + l)\d0. 
Jo 

(18) 

For the case of primary interest to us, m = 2, this integral 
has the value X 2 = 0.80411 [4]. 

Now we can define the viscous timescale for dissipation 
in the boundary layer: 



1 



1 dE 
2E dt ' 



(19) 



Using the expression for dE/dt derived above, and the 
usual expression for the energy E we find: 



1 2 m+3 / 2 (m + l)! 
Tv ~ 20m(2ra+l)!!I ra X 



p_ 

Pc 



2m+2 



dr 

Rc 



(20) 



Here the quantities R c , p c and r\ c are the radius, density, 
and the viscosity of the fluid at the outer edge of the 
core. We note that, while the viscous dissipation rate in 
the boundary layer, Eq. (17), depends only on quantities 
local to the boundary layer (d, p c , R Cl CI), the timescale 
t v depends also on the global structure of the mode. This 
is due to the fact that, while most of the energy dissipa- 
tion takes place in the boundary layer, most of the energy 
in the mode is not localized there. For realistic neutron 
stars this expression, Eq. (20), for the viscous timescale 
is about a factor of two larger than the one obtained 
by Rieutord [4], who assumed a uniform-density stellar 
model. 

In deriving our expression for the viscous boundary 
layer dissipation timescale, Eq. (20), we assumed that the 
crust is rigid and hence static in the rotating frame. The 
motion of the crust due to the mechanical coupling to the 
core [6] effectively increases t v by a factor of (Av/v)~ 2 . 



IV. STABILITY OF THE fl-MODES 

In this section we evaluate viscous timescales for neu- 
tron stars (based on slowly rotating Newtonian models) 
constructed from a set of "realistic" equations of state, 
and use these timescales to evaluate the stability of these 
stars. In neutron stars colder than about 10 9 K the shear 
viscosity is expected to be dominated by electron-electron 
scattering. The viscosity associated with this process is 
given by [8,9] 



rf e = 6.0 x 10Vr 



(21) 



where all quantities are given in cgs units, and T is 
measured in K. For temperatures above about 10 9 K, 
neutron-neutron scattering provides the dominant dissi- 
pation mechanism. In this range the viscosity is given by 
[8,9] 



r,™" = M7p 9/4 T- 



(22) 



We find it useful to factor the angular velocity and tem- 
perature dependence from the viscous timescale defined 
in Eq. (20). Thus we define a fiducial viscous timescale 
f v such that 



n 



1/2 ( T \ 



(23) 



where VL Q = \J-nGp. We have evaluated these fiducial vis- 
cous timescales (for each type of viscous dissipation) for 
1.4M neutron star models based on a variety of realis- 
tic equations of state [10] as well as the standard n = 1 
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TABLE I. Properties of 1.4M neutron stars with rigid 
crusts for densities below p c — 1.5 x 10 14 g/cm 3 . Times are 
given in seconds, lengths in kilometers, temperatures in units 
of 10 8 K, and energies in units of 10 51 ergs. 



EOS a 



R 



TGR 



rpj} 



nlpoly 
BJI 
DiazII 
FP 

Glendl 

Glend2 

HKP 

PandN 

SHW 

WFF3 

WGW 



12.53 
15.04 
15.80 
13.18 
16.47 
16.84 
16.43 
12.77 
14.92 
12.98 
14.64 



11.01 
12.16 
12.60 
11.34 
12.74 
12.95 
13.27 
10.85 
12.89 
11.21 
12.38 



1.94 
1.30 
1.20 
1.93 
0.97 
0.96 
1.32 
1.84 
1.97 
2.10 
1.85 



2.53 
1.72 
1.67 
2.18 
1.53 
1.53 
1.76 
2.08 
2.23 
2.33 
2.17 



4.25 
13.11 
17.37 
5.24 
25.32 
28.12 
18.41 
4.84 
8.40 
4.52 
8.32 



23.3 
24.6 
24.3 
25.9 
23.4 
23.7 
24.6 
26.9 
25.0 
27.7 
26.8 



1.69 
4.96 
6.64 
1.88 
10.1 
11.0 
6.97 
1.67 
3.11 
1.52 
2.88 



52.0 
54.8 
54.2 
57.8 
52.2 
52.8 
54.8 
59.9 
55.8 
61.8 
59.8 



0.76 
2.23 
2.98 
0.84 
4.51 
4.95 
3.13 
0.75 
1.40 
0.68 
1.29 



a The various equations of state used here are described in 
Bonazzola, Frieben and Gourgoulhon [10]. 
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FIG. 1. Temperature dependence of the critical angular ve- 
locities of 1.4M Q neutron stars constructed with a number of 
realistic equations of state. 



polytrope with a radius of 12.53 km. These results are 
summarized in Table I, along with other relevant proper- 
ties of these stellar models. In particular we also include 
the total radius R, the radius of the fluid core R Cl the en- 
ergy contained in an excited r-mode e = E(Q /aQ) 2 ergs 
for the case of a fully fluid star, the energy contained in 
an excited r-modc e c in the case with rigid crust, the 
fiducial gravitational radiation timcscale tgr, and a crit- 
ical temperature T c (defined below). We note that the 
t v values listed in Table I presume that the crust is rigid 
and does not move in the corotating frame. To take into 
account the motion of the crust, multiply f^ e and f™ by 
(Av/v)- 2 [11]. 

The gravitational radiation timescale is evaluated here 
using the formalism developed by Lindblom, Owen, and 
Morsink [12]. This timescale is given by 



1 _ 327rGfi 2m + 2 (to - l) 2m 

TGR 



c 2m+3 



X 



TO. - 



[(2to + 1)!!] 2 

2m+2 r R c 



1 



pr 



2m+2 



dr. 



(24) 



Since we presume that the crust is static in the rotat- 
ing frame, the integral in Eq. (24) extends only over 
the interior of the fluid core, rather than throughout 
the star. Consequently, these gravitational radiation 
timescales are somewhat longer (typically about 30% 
longer) than those computed earlier for hotter, com- 
pletely fluid stars [11]. The fiducial gravitational radi- 
ation timescale t G r given in Table I is defined by 



vl \ f v J 



2/H / 10 8K X 2/11 



T 



(26) 



Figure 1 illustrates the temperature dependence of the 
critical angular velocity for 1.4M Q neutron stars con- 
structed from a variety of realistic equations of state. 
The discontinuities in these curves at T = 10 9 K occur 
because the supcrfluid transition changes the viscosity 
from electron-electron scattering at low temperatures to 
neutron-neutron scattering. The outer edge of the core 
in these models is taken to be p c = 1.5 x 10 14 g cm~ 3 [13]. 
Equation (26) and the curves in Fig. 1 neglect the mo- 
tion of the crust in the rotating frame, i.e., they assume 
that Av/v = 1. In general, O c is a factor of {Av/v) 4 / 11 
smaller than that given in Eq. (26). Thus, fl c displayed 
in Fig. 1 provides an upper bound on the critical fre- 
quency for realistic neutron stars. In these expressions 
for fl c we neglect the effects of bulk viscosity, which are 
unimportant for T < 10 10 K. 

The angular velocity of a neutron star can never exceed 
some maximum value f2 max sa |f2 Q [10]. Thus, there is 
a critical temperature below which the gravitational ra- 
diation instability is completely suppressed by viscosity. 
This critical temperature is given by 



10 8 K 



11/2 - 

TGR 



gX 11/2 ~ 

° 1 ^P. (27) 



In terms of T c then the critical angular velocity can be 
expressed in a particularly simple form: 



tgr = tgr 



2m+2 



(25) 



n 



fin 



2/11 



2/11 



(28) 



Gravitational radiation tends to drive the r-modes un- 
stable, while viscosity suppresses the instability. We de- 
fine the critical angular velocity f2 c , above which the r- 
mode is unstable, by the condition t v = tgr: 



The values of this critical temperature are given in Ta- 
ble I for the case p c = 1.5 x 10 14 g cm' 3 . Since the exact 
density p c where crust formation begins is only poorly 
known, we explore in Fig. 2 the dependence of T c on p c . 
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p c /10 14 g/cm 3 

FIG. 2. Critical temperature T c as a function of the crust 
formation density p c for a number of realistic equations of 
state. 

The motion of the crust would change T c of Eq. (27) by 
the factor (Av/v) 2 , and so the curves in Fig. 2 provide 
an upper bound on T c for realistic neutron star models. 



V. THERMAL STRUCTURE OF THE 
BOUNDARY LAYER 

Viscous dissipation in the boundary layer between the 
core and the crust deposits thermal energy into a very 
thin layer of material. Previous authors assumed that 
this heat is effectively conducted away from the bound- 
ary layer, so that the star remains essentially isothermal. 
Clearly, this situation is idealized, and viscous dissipa- 
tion will raise the temperature at the crust-core inter- 
face. The viscous heating competes with thermal con- 
duction away from the boundary layer, as well as neu- 
trino emission from the crust and the core. In this sec- 
tion we evaluate the r-mode amplitude needed to raise 
the temperature at the boundary layer to a given value. 
Our particular interest is to determine the minimum 
r-mode amplitude needed to raise the temperature to 
T m w 10 10 (/9/p c ) 1 / 3 K, the melting temperature of the 
crust (see the Appendix for a discussion of the micro- 
physics employed in this paper). If this occurs, the r- 
modes may continue to be unstable to much lower angu- 
lar velocities (see Sec. VI). 

In the discussion that follows, it will be necessary to 
know the explicit temperature dependences of the ther- 
mal conductivity and neutrino emissivity. As described 
in the Appendix, these temperature dependences are 



k = 



~rpr, 



(29) 



(30) 



where T\q = T/10 10 K, and the exponent in the emissiv- 
ity equation is n = 8 in the fluid core and n = 6 in the 



crust. The prefactors k and e are independent of temper- 
ature and are given by k< = 1.5 x 10 31 , e< = 8.6 x 10 28 
in the core and k> = 2.8 x 10 30 , e> = 1.5 x 10 25 in the 
crust. The quantities n and e are in cgs units. Note that 
for notational convenience in subsequent computations, 
the temperature in Eq. (29) is measured in K, while in 
Eq. (30) it is measured in 10 10 K. 

Let us first neglect the thermal conductivity alto- 
gether, and ask what r-mode amplitude is necessary to 
heat the boundary layer to T m if the heat is radiated ex- 
clusively by neutrino emission from the boundary layer 
itself. The rate dE/dt at which the shear deposits en- 
ergy into the vicinity of the boundary layer is given by 
Eq. (17), where the amplitude A is related to the di- 
mensionless r-mode amplitude, as defined by Lindblom, 
Owen and Morsink [12], via 



a 



IQttR 



A. 



(31) 



We equate this rate to the neutrino emission rate in the 
boundary layer, A-k i? 2 iie < T 1 8 , to obtain a c (local), the r- 
mode amplitude necessary to keep the crust-core inter- 
face at the melting temperature if thermal conduction is 
not important: 



Qc (local)=-^ 5W20 3 j 



1/2 



(32) 



In this equation the numerical prefactor has been evalu- 
ated using the standard n = 1 polytropic stellar model. 
This estimate sets the lower bound on the critical melting 
amplitude, since heat conduction increases the volume 
that radiates heat by neutrino emission, and hence more 
viscous dissipation is required to maintain the interface 
at a certain temperature. 

We now show how to evaluate the effect of finite con- 
ductivity. The general equation for the thermal evolution 
of the material in a neutron star is 



C v 8 t T = -e + V • (kVT) + 2 V Sa* ab a 



ab 



(33) 



where T is the temperature, C v is the specific heat at con- 
stant density, e is the neutrino emissivity, and k is the 
thermal conductivity. This is a time-dependent equation, 
and in general the overall cooling of the star due to Urea 
neutrino emission must be followed along with the heat- 
ing in the boundary layer. Let £ denote the thickness of 
the region adjacent to the boundary layer whose temper- 
ature is raised above the ambient by the viscous dissipa- 
tion. On a timescale idiff ~ C v £ 2 /k the heat generated in 
the boundary layer will diffuse throughout this larger re- 
gion. In the neutron-star matter near the boundary layer 
this timescale is idiff ~ 8 x 10~ 5 (£/d) 2 s, which depends 
on the temperature only through the boundary layer 



G 



thickness, d oc T -1 . The neutron star as a whole cools on 
the Urea cooling timescale, t coo \ ss 30 T^ 6 s [14]. If the 
diffusion timescale is shorter than the cooling timescale, 
then the temperature distribution in the region near the 
boundary layer will estabilish a quasi-equilibrium state in 
which the excess heat is conducted into a large enough 
volume for it to be radiated away by neutrinos. The tem- 
perature at the inner edge of this region slowly decreases, 
tracking the overall cooling of the star. As we shall see 
below, the width of this quasi-equilibrium layer is a few 
hundred times the thickness of the boundary layer itself. 
Thus idiff < ^cooi and so a quasi-equilibrium state exists 
in which C v dtT can be neglected compared to the heat 
conduction term d r (nd r T). 

In the quasi-equilibrium region, but outside the bound- 
ary layer, the temperature distribution is described ap- 
proximately by 



d r {Kd r T) = e. 



(34) 



Using the simple dependence of n and e on temperature, 
given by Eqs. (29) and (30), it is straightforward to ob- 
tain a first integral of Eq. (34), 



dT w 
dr 



2e 



UK 



rp2 frpn 

~ 1 io v J 10 



T, 



o.lOJ' 



(35) 



where T Q is the ambient temperature in the core where 
the heat flux nd r T tends to zero. 

The viscous energy generation in the boundary layer 
is determined by the shear given by Eq. (14). Here we 
re-express this energy generation rate as 



F(r,0), 



(36) 



where the function F(r,0) is defined in Eq. (15). This 
function falls off exponentially away from the boundary 
layer with the length scale d. Thus F can be reasonably 
approximated as a delta function, 



F(r,6) 



477 S(r-R c ) 



sn Pc p{6) ' 

where the angular function f(6) is defined by 



(37) 



f- 2 (6) = J % sin 2 6 (1 - cos6») 2 p+ + (l + cos0) 2 p_ 



(38) 

Now we return to the full equation for the thermal 
distribution in the vicinity of the boundary layer: 



d r {nd r T) =e-2r]Sa* b Sa 



ah 



(39) 



This equation can be integrated analytically, using 
Eq. (35), when we approximate the heating in the vis- 
cous boundary layer by Eq. (37) . The r-mode amplitude 
a c needed to raise the temperature to the value T m is 
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FIG. 3. Dependence of the critical amplitude a on T a for 
stars with fl = 0.5fi o . The smooth curve is the analytical 
expression given in Eq. (40), while the points are numerical 
solutions to the heat equation without the delta-function ap- 
proximation for the boundary layer heating term. 



2 bR^ijp^ _ pr — - 
c 16VZnRip(6) ~ V 4 £ < K <i J m ,io *o,io) 

V|e>«>(T6 il0 -r o 6 il0 ). (40) 



For the values of the microphysical parameters described 
above, the critical amplitude satisfies 



^,o vd =8 - 0xl ° T"5t 



l 5-3 x 10-8 T6 

+ ~f^rf'n l - (41) 



We illustrate in Fig. 3 the dependence of this critical 
amplitude on the ambient temperature T Q for the case 
T m = 10 10 K and 8 = tt/2. The solid curve is the an- 
alytical expression given in Eq. (41). For comparison, 
we also include a selection of points computed numeri- 
cally by solving the full differential equation, Eq. (39), 
without making the assumption that the heat source is a 
delta function. Clearly, the analytical approximation is 
extremely good. We also see that, because of the steep 
dependence of the neutrino cooling rate on the temper- 
ature, the value of a c is rather insensitive to T a , except 
when T fa T m . 

The prefactor of the second term on the right side of 
Eq. (41) is much smaller than the prefactor of the first 
term. Thus we see that the core plays a much more im- 
portant role in determining a c than the crust. This is 
easy to understand: At T = T m , the neutrino emissiv- 
ity in the core is approximately 5000 times larger than 
in the crust. Since we demand that the crust-core in- 
terface is kept at a fixed temperature, much more heat 
is emitted on the core side of the interface than on the 
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cose 

FIG. 4. Function / that determines the angular depen- 
dence of the critical amplitude a. 

crust side. Therefore, the heat flux into the core must be 
much higher than the flux into the crust. Since the heat 
flux into the crust can be neglected (i.e. e<«< 3> e>K>), 
we may obtain a much simpler expression for the critical 
amplitude: 

a c » 2.8 x lO-^S/W • (42) 

for the case T a -C T m . For a maximally rotating neu- 
tron star, with f2 max w |0 o , the prefactor in the above 
equation is 4.7 x 10~ 3 . 

The function f(0), illustrated in Fig. 4, determines 
the angular dependence of the critical melting amplitude. 
Melting the crust near the poles requires a higher r-mode 
amplitude than melting at the equator. However, over a 
wide range of angles the critical amplitude needed to heat 
the interface to a given temperature is the same as the 
equatorial (9 = w/2) value to within a factor of two. 

The critical melting amplitude a c is roughly 20 times 
a c (local), our estimate that neglects the effect of ther- 
mal conduction, Eq. (32). From this we may deduce 
that the thickness of the layer into which thermal en- 
ergy is conducted before being radiated away is of order 
II dpi a 2 c /a 2 (local) w 400. From our numerical solutions, 
we confirm that the thermal flux in the core falls from 
its value at the crust-core interface to half that value at 
approximately 500c?. Thus the thickness of the layer that 
radiates the excess thermal energy is small enough (just a 
few centimeters) to justify the quasi-equilibrium analysis 
which we used. 

Finally, we note that the heating rate used in this cal- 
culation neglects the motion of the crust in the rotating 
frame [6] . The motion of the crust increases the critical 
melting amplitude by the factor of (Av/v)^ 1 . 



VI. DISSIPATION IN THE ICE FLOW 

As a young neutron star cools, its temperature quickly 
falls below the freezing temperature and a crust begins 
to form. If no unstable r-mode is present, a solid crust 
forms as usual. However if the star is rapidly rotating, 
the m = 2 r-mode will be driven unstable, and may grow 
to an appreciable amplitude even before the crust begins 
to form. Therefore, the crust would have to form in the 
presence of the shearing motion of the r-mode. Even in 
the presence of the crust, the r-mode is unstable at high 
enough frequencies, as illustrated in Fig. 1. Once the am- 
plitude of the mode grows beyond the critical value given 
in Eq. (42), the dissipation in the boundary layer would 
be sufficient to re-melt the crust if it formed. If the crust 
fails to form or is completely melted away, the boundary 
layer heating would dissappear and the dissipation would 
not be sufficient to prevent freezing. Clearly, the state of 
the outer layers of the star can be neither pure solid nor 
pure fluid. We imagine that instead a mixed state will 
form: a neutron-star ice flow, much like the pack ice on 
the surface of the Arctic Ocean. 

Microscopic viscosity will be replaced in this ice flow 
by dissipation that is dominated by collisions between 
macroscopic chunks of crust, the boundary layers around 
the ice chunks, or other mechanisms. However, regard- 
less of the detailed dissipation mechanism, this flow is 
self-regulating: as long as the amplitude of the r-mode 
remains above the critical value of Eq. (42), the dissipa- 
tion must, on average, be enough to keep the crust at 
around the melting temperature. Any less dissipation, 
and the crust would completely freeze, thus giving rise 
to a boundary layer which would re- melt the crust. Any 
more dissipation and the crust would simply melt, then 
promptly re-form as the fluid viscosity alone is insuffi- 
cient to keep the outer layers hot. The assumption of 
self-regulating flow allows us to compute the dissipation 
rate regardless of the details of the dissipation mecha- 
nism. 

The dominant energy sink at high temperatures is 
bremsstrahlung neutrino emission, with cmissivity given 
by Eq. (A6). Thus the viscous dissipation in the ice flow 
must just balance the losses due to neutrino emission. 
Presuming that the outer layers of the neutron star are 
kept at the local melting temperature given by Eq. (A2), 
we can determine the total energy dissipated in the ice 
flow, and hence an effective timescale Ti ce on which this 
ice flow dissipation will damp the r-mode: 

— = ^[eT^ w d 3 x. (43) 

Displaying explicitly the scalings of t; co with the mode 
amplitude a and the angular velocity fi, we find 

Ticc = a 2 (JpJ r ico . (44) 
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The fiducial timescales fi ce range from 1.4 x 10 8 s to 
1.3 x 10 9 s for the "realistic" equations of state discussed 
above. Unlike the other dissipative timescales in this 
problem, this one depends on the mode amplitude a. 

The r-mode will continue to be unstable, and thus will 
continue to spin down the neutron star by emitting grav- 
itational radiation, as long as the gravitational radiation 
timescale is shorter than the ice flow viscosity timescale, 
and its amplitude remains above the value needed to sus- 
tain the ice flow. The critical angular velocity where the 
equality is achieved and the mode becomes stable is 



1/8 



o (—) ^ ' 



(45) 



For a = 1, this critical angular velocity has the value 
0.093O o for the standard n = 1 polytropic stellar model, 
and values that range from 0.086f2 o to 0.14O o for the real- 
istic models considered here. (Note that the gravitational 
radiation timescale tqr used here is the one appropriate 
for a mode that extends to the surface of the star, not 
the one given in Table I for a mode that extends only to 
the edge of the core.) Thus we see that if the amplitude 
of the r-mode saturates at a value that is close to unity, 
the instability will continue — even in the presence of a 
crust — until the angular velocity of the star is a small 
fraction of il . 

Is the ice-flow picture described above consistent? 
Clearly the dissipation due to molecular viscosity is in- 
sufficient to keep the outer layers of the neutron star 
hot, but can collisions between macroscopic ice chunks 
accomplish this? Let us estimate the characteristic size 
of the ice chunks needed to produce the needed macro- 
scopic viscosity. We assume that the ice chunks are of 
average size D, and the mean free path between them 
is A. If A D, then the collisions between the chunks 
are rare, and do not significantly raise the viscosity com- 
pared to the purely fluid value (see § 22 of Landau and 
Lifshitz [15]). We therefore consider the regime where 
the ice chunks are densely packed, with D > A, and be- 
have like a fluid. We follow the discussion of Haff [16] 
and Borderies et al. [17]. 

The ice chunks mainly follow the r-mode velocity flow, 
but due to the collisions between themselves acquire a 
random component of velocity, v. The viscosity in such 
a flow is given by ?7i cc « pD 2 v/X [16,17], and is a factor 
of D/X bigger than the usual molecular viscosity. This 
is because each particle transports momentum across a 
distance D while only having traveled the (possibly much 
smaller) distance A between each collision. 

The collisions between the chunks are inelastic with 
a fraction 7 of the collision energy per chunk, pD 3 v 2 , 
being converted to heat. Considerations of icy particle 
collisions in the temperature regime appropriate for Sat- 
urn's rings put 7 « 0.7 — 0.9 (see Borderies et al. [18] for a 
review). Since our flow is near the melting temperature, 
we expect very inelastic collisions as well, i.e., 7 w 1. 
The collision rate is just v/X, so the energy dissipation 



per unit volume in the flow is jpv 3 /X. For the r-modes 
the shear dissipation rate is 2r]i ce 5al b 5<T ab w r]i cc a 2 fl 2 , 
so v w 7 _1 / 2 aOD. Thus when the collisions are inelas- 
tic, the random component of the ice chunk velocity is 
comparable to the velocity difference between the neigh- 
boring chunks in the overall ice flow. The viscosity of 
this granular flow is therefore given by 



'/i< 



paQD 3 
1 / 2 X 



10 



18 



a 



D 3 



7 



1/2 A 



7 



n 
n~ 



(46) 



where A and D are measured in centimeters and 77i ce 
is in cgs units. Not surprisingly, this viscosity is much 
larger than the microscopic viscosities given by Eqs. (21) 
and (22). 

The energy dissipation rate in the ice flow is obtained 
by integrating rii cc 8(T* ab 8(T ab over the outer layers of the 
star: 



dE 
~dt 



(attD) 3 M c 



7 



(47) 



where M cr is the mass of the material at densities less 
than p c . Using the r-mode energy for an n = 1 poly- 
trope (see Table I), the damping time due to this form 
of viscous dissipation is 



10 8 s 



7 



1/2 A 



aD 3 



0.05 1/ j 



(48) 



The heat deposited by the ice flow into the star is radi- 
ated by neutrinos on the timescale evaluated in Eq. (44). 
Equating these two timescales gives an estimate of the 
average ice chunk size necessary to keep the crust at the 
local melting temperature: 



10 8 s 



1/2 A 



D 



0.05MA 



1 1/2 



(49) 



Our estimate of the ice chunk size depends only weakly 
on the unknown inelasticity of the collisions and the ratio 
X/D. Wc do not expect A to be much smaller than D, 
because if it were, the ice chunks would probably lock 
together into bigger pieces, leading to increased friction 
and dissipation. Moreover, as argued above, we expect 
7 w 1 as well. Thus for typical values of fi ce evaluated 
above, the value of the ice chunk size needed to keep 
the outer layers of the star at the melting temperature is 
D w 1 cm. 

Several conclusions can be drawn from this estimate. 
First, the smaller the mode amplitude, the larger the ice 
chunks have to be in order to provide enough friction 
(recall that the viscosity is proportional to D 2 ) to keep 
the fluid at the melting temperature in the face of neu- 
trino losses. The chunk size exceeds the radius of the star 
when a < 10- 4 (n o /Q)(X/ D) 1 / 3 . Hence, if the mode am- 
plitude is smaller than this value, the ice-flow mechanism 
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is unable to keep the outer part of the star in a melted 
state, and it must freeze. Since this amplitude is smaller 
than the critical melting amplitude computed in the pre- 
vious section, this does not further restrict the range of 
mode amplitudes where the ice flow is expected to occur. 
Second, even for a w 1 the chunk size is small but not 
microscopic. So our picture of the ice flow is consistent 
throughout the regime of interest. 

The argument that the ice flow maintains the star's 
outer layers near melting depends on the ability of the r- 
mode to re- heat a solid crust back to the verge of melting, 
should a crust form through some fluctuation. Up to this 
point we assumed that this heating is instantaneous once 
a exceeds a c . This is not an issue if f2 > fl c [the criti- 
cal angular velocity for the r-mode instability, Eq. (26)], 
since the amplitude of the r-mode is growing due to ra- 
diation reaction in this case. However, if Q < f2 c , the 
r-mode is decaying due to the dissipation in the bound- 
ary layer even while it is re- heating the crust. Can the 
r-mode re-heat the crust and re-form the ice flow be- 
fore being damped? A simple comparison of the energy 
needed to melt the crust, E m , to the mode energy, E (see 
Table I), is inadequate, since it neglects the continuing 
energy input from gravitational radiation reaction. 

The r-mode melts the crust on a timescale 



E„ 



En 



7~rn — 



(3(dE/dt) v 2f3E 



(50) 



where the factor (3 is the fraction of heat flowing 
into the crust, which is given approximately by (3 = 
(4? > K > /3e < K < ) 1 / 2 w 0.007 (see Eq. [35]). While it is 
melting the crust, the r-mode is also damped at the rate 
1/Td = 1/t v — 1/tqr. The damping rate is exactly zero 
on the stability curve = fl c . However, since the r-mode 
has a finite amplitude, it will continue to spin the star 
down below the stability line while melting the crust. 

Let Afi be the change in spin frequency, due to the 
r-mode evolution, during the time r TO . For small a, the 
spindown rate is given by Eq. (3.14) of Owen et al. [19]: 



1 dn 

Q~dt 



0.2a 2 



(51) 



Thus, in the time it takes to melt the crust, the star spins 
down by Af2/Q c w 0.1a 2 E m / f3E (which is independent 
of a since E ca a 2 ). By using a Taylor expansion, it 
is easy to show that the r-mode damping time for Q, = 
Q c — Afl is approximately 



2 n c 
Td ~iiMi Tv 



2/3E 
a 2 E, 



-Ty. 



(52) 



If T m < Td, we conclude that the r-mode can re- melt the 
crust and so re-create the ice flow before being damped 
out. This condition places the following limit on the 
mode amplitude a, 



a 



> 



Em 

2f3e 



3 x 10" 3 



(53) 



In the above equation, we assumed that the energy re- 
quired to re-melt the crust is of order kT per nucleus 
(since the crust is near the melting temperature), i.e. 
E m w 10 47 T m ,io erg. This amplitude is very close, both 
in the absolute magnitude and in the scaling with and 
T m , to the critical melting amplitude a c . Hence, so long 
as the r-mode amplitude is large enough to raise the 
crust-core interface to the melting temperature, the en- 
ergy contained in the mode is enough to melt the crust 
even if the r-mode is no longer linearly unstable. Since 
the criterion (53) does not depend on Av/v while a c 
(Eq. [42]) increases with Av/v, the conclusion does not 
change even when the motion of the crust is taken into 
account. 



VII. DISCUSSION 

We have computed stability curves (critical frequency 
as a function of temperature) for the onset of the r- 
mode instability in a neutron star with a (laminar) vis- 
cous boundary layer under a solid crust. We improve 
previous calculations by including the effect of the Cori- 
olis force (the dominant restoring force) on the bound- 
ary layer and by using realistic neutron-star models — two 
important ingredients which have not been combined in 
previous work. Our stability results are summarized in 
Eqs. (27) and (28), Table I, and Figs. 1 and 2. If the 
neutron star crust is rigid and does not move in the ro- 
tating frame, then our results imply that the r-modes 
are not unstable in any of the accreting neutron stars in 
LMXBs or millisecond pulsars. However, if the relative 
velocity amplitude between the core and the crust, Av/v 
is significantly smaller than 1, as recent calculations for 
constant-density stars suggest [6], then our results con- 
stitute only the upper limits on the critical frequencies 
and temperatures for the onset of the r-mode instability. 
A self-consistent calculation of the r-mode eigenfunctions 
in the presence of a realistic crust is necessary to settle 
(at least for unmagnetized neutron stars) the question of 
linear stability of the r-modes. 

We find that localized heating in the boundary layer 
between the solid crust and the fluid core can successfully 
compete both with the heat conduction away from the 
boundary layer and with neutrino emission. In Section V, 
we computed the critical r-mode amplitude, Eq. (42), 
needed to raise the temperature of the boundary layer 
to the crust melting temperature while the interior of 
the star far away from the boundary layer remains at a 
much lower temperature. The amplitude required for this 
is rather small, of order 10~ 2 , and is comparable to that 
necessary to break up the crust via purely mechanical 
coupling. 

Based on the smallness of the melting amplitude, we 
argue in Section VI that the spindown scenarios for neu- 
tron stars with crusts advanced by previous authors may 
need to be significantly modified. So long as the r-mode 
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amplitude remains above the melting value, a completely 
solid crust is not possible. Instead, the outer layers of a 
neutron star will resemble an ice pack which flows along 
with the r-mode motion. Viscosity in the ice pack ad- 
justs itself to maintain the outer layers of the neutron 
star near the crust's melting temperature, regardless of 
the details of the dissipation in the ice flow. The ice flow 
can persist even if the r-mode is no longer linearly un- 
stable according to the instability criterion for neutron 
stars with crusts, since an r-mode with an amplitude 
greater than the minimum value given in Eq. (53) can 
melt the crust before being damped by boundary-layer 
friction. Therefore, if the crust is melted cither during 
spindown of young neutron stars, or in the final stage of 
a thermal runaway in LMXBs, the final spin frequency 
is not set by the boundary-layer damping time, but in- 
stead by the balance between n cc and tqr. This leads 
to much smaller final spin frequencies than computed by 
previous authors. A neutron star with an r-mode ampli- 
tude large enough to create the ice pack can spin down 
to a frequency nearly as low as a purely fluid star [see 
Eq. (45)], provided the saturation amplitude of a purely 
fluid r-mode of is order unity. 

Very recently, Wu et al. [21] evaluated the effect of 
a turbulent boundary layer between the crust and the 
core on the stability of the r-modes. They find that 
the turbulent dissipation timescale depends on the r- 
mode amplitude and therefore determines a saturation 
amplitude, which they estimate (in our notation) as 
a sat w 1.5 x 10- 2 (n/n o ) 5 (Av/v)- 3 , or about 2.0 x 10~ 3 
for a maximally rotating star with a rigid crust. Inserting 
Eq. (10) of Wu et al. [21] into our Eq. (39), we find that 
a turbulent boundary layer changes the critical melting 
amplitude at the equator to 



(54) 



or a c 7.9 x 10~ 4 for a maximally rotating polytropic 
neutron star with a rigid crust. Comparing with Eq. (42), 
we conclude that the crust can be heated to the melting 
temperature at even smaller amplitudes than in the pres- 
ence of a laminar boundary layer. Thus, if the r-mode 
amplitude in a newborn neutron star exceeds the critical 
value given above when the temperature drops below the 
crust's melting temperature, then crust formation will be 
delayed (as we argued in Sec. VI). In this case the star 
will spin down to frequencies close to those predicted for 
crustless fluid stars, provided the saturation amplitude 
of the r-mode is of order unity. 
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APPENDIX: THERMODYNAMIC PROPERTIES 
OF NEUTRON-STAR MATTER 

The strength of electrostatic interactions between the 
nuclei in neutron-star matter is typically expressed in 
terms of the Coulomb coupling parameter, 



2 2 



Z 2 e 
akT ' 



(Al) 



where Z is the nuclear charge, and a = (3A ccll m;,/47rp) 1 / 3 
is the radius of the Wigner-Seitz cell with A ce n nucle- 
ons in it. According to the recent molecular-dynamics 
simulations of Farouki and Hamaguchi [22] , such a clas- 
sical one-component plasma crystallizes when T > 173. 
However, the details of the nucleon-nucleon interaction 
that determine Z and A ce n near the bottom of the crust, 
p = p c , are still not completely understood. Depending 
on the particular interaction model, values of Z ranging 
from approximately 30 to 50, and A cc \\ in the range of 
500 to 1000 are obtained. These values result in melt- 
ing temperatures within a factor of two of 10 10 K. (Sec 
Pethick and Ravenhall [23] for a review of the nuclear 
physics, Douchin and Hacnsel [24] for the latest cal- 
culation, and Haensel [25] for a discussion of melting 
temperatures appropriate for young neutron stars and 
LMXBs.) Moreover, for some nuclear force models, nu- 
clei at p > 10 14 g cm~ 3 may resemble rods, plates, and 
tubes, rather than spheres (Lorentz et al. [26]). Melting 
temperatures of such exotic phases of matter have not 
been calculated, but it is reasonable to assume [27] that 
they will be comparable to those of "ordinary" matter 
with spherical nuclei. In the light of these uncertainties, 
we adopt 



T m = 10 



10 



1/3 



K, 



(A2) 



as the fiducial melting temperature of crustal matter. 

The thermal conductivity in the neutron-star core is 
dominated by electron-electron collisions. A convenient 
fit to the conductivity is given by Flowers and Itoh [28] : 



1.5 x 10 21 



erg cm 



1 K" 



(A3) 



at p = p c , where T w = T/10 10 K. The neutrino emis- 
sivity in the core is given by the usual modified Urea 
formula [29,30]: 
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e = 8.6 x 10 28 T 8 erg cm- 3 s" 1 , (A4) 

at p = p c . These lead to the expressions for the core 
values, k< = 1.5 x 10 31 and e< = 8.6 x 10 28 , used in 
Sec. V. 

The microphysics in the crust is, in general, quite a bit 
more complicated. However, we are primarily interested 
in the conditions at or near the melting temperature. For 
T > 2 x 10 10 K, the crustal conductivity is also dominated 
by electron-electron collisions, and even at T 10 10 K, 
to within a factor of two, 

2 - 8 X 1020 -1 -1 T^-l { \ K\ 

k w n ee w — erg cm s K , (A5) 

J-w 

where we approximated the fits of Flowers and Itoh [28] 
for the temperature regime of interest. For temper- 
atures lower than this, the conductivity is dominated 
by electron-phonon scattering, which is approximately 
constant, « c -ph ~ 10 20 for temperatures in excess of 
the Debye temperature, T d w 5 x 10 9 (p/p c ) 1 / 2 K. For 
T <C Td, the electron-phonon scattering has the same 
temperature dependence (but a different prefactor) as 
electron-electron scattering, K c _ p h ~ 1 x 10 19 /Tio- We 
also find that the recent calculations of neutrino-pair 
bremsstrahlung in the crust (Hacnscl et al. [31] and 
Kaminkcr et al. [32]) can be reasonably approximated 
by 

e w 1.5 x 10 25 (j^J T\ 6 erg cm" 3 s" 1 . (A6) 

These expressions lead to the crust values, k> = 2.8 x 
10 30 and e> = 1.5 x 10 25 that we use in Sec. V. Since the 
crust does not play a significant role in determining the 
temperature at the boundary layer (as shown in Sec. V), 
including a detailed treatment of the microphysics there 
is not essential. 
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